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Abstract 

A search for Higgs bosons produced in association with a fermion pair, and 
decaying to WW, is performed with the data collected by the ALEPH detector at 
centre-of-mass energies ranging from 191 to 209 GeV. The data correspond to an 
integrated luminosity of 453.2 pb _1 . Thirteen exclusive selections are developed 
according to the different final state topologies. No statistically significant evidence 
for a Higgs boson decaying into a WW pair has been found. An upper limit is 
derived, as a function of the Higgs boson mass, on the product of the e + e~ — ► 
Hff cross section and the H —* WW branching ratio. The data on the search for 
H — > WW are combined with previously published ALEPH results on the search 
for H — > 77, to significantly extend the limits on the mass of a fermiophobic Higgs 
boson. 
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1 Introduction 



A Higgs model pQ incorporating two doublets of complex scalar fields [2j generates five 
scalar Higgs bosons, three of which are neutral. In some types of models, for certain choices 
of parameters, one of these neutral scalars provides mass only to the fermions and another 
couples exclusively to the bosons, i.e. is a "fermiophobic" Higgs boson. Anomalous 
couplings in the Higgs sector can also enhance the bosonic branching fraction 

The search for a fermiophobic Higgs boson has been primarily carried out by the four 
LEP experiments in the H — > 77 channel, in which the Higgs boson couples to photons 
via a W loop [U El El [7]. A benchmark fermiophobic Higgs boson is defined by considering 
Standard-Model-like couplings to bosons, and null couplings to fermions. Current analyses 
exclude the benchmark fermiophobic Higgs boson up to a mass of 109.7 GeV/c 2 [8j. For 
fermiophobic Higgs bosons heavier than 90 GeV/c 2 , the predicted H — > 77 branching ratio 
becomes small relative to the predicted H — > WW branching ratio (Fig. motivating a 
search in this new channel. Such an analysis has already been carried out by the L3 
collaboration JU] and is performed here with data collected by the ALEPH detector. 

The main production processes at e + e~ colliders for a fermiophobic Higgs boson are 
e + e~ — > Z* — > ZH (Higgsstrahlung), WW and ZZ fusion. The cross sections of the boson 
fusion production processes are considerably smaller than that of the Higgsstrahlung 
process at LEP centre-of-mass (CM) energies. In the mass range kinematically accessible 
for Higgsstrahlung at LEP, one of the virtual W bosons is expected to be near on-shell, 
and the other (denoted W*) to have a much smaller mass and energy. In this paper, all the 
signatures originating from the Z — > qq, z/z/, £ + £~ decays and the W, W* — > qq', £ ± z/ decays 
are searched for. For simplicity and conciseness, the term "lepton" (and the corresponding 
symbol £) refers to electrons and muons only. Leptonic tau decays are not specifically 
addressed, but the corresponding selected events are included in the final results. The 
hadronic tau decays in e + e~ — > ZWW* — > £ + £rrvc^[ are also looked for. The analysis is 
performed on the data taken in the years 1999 and 2000 at CM energies ranging from 191 
to 209 GeV. The luminosities and CM energies are shown in Table H 

This paper is organized as follows. A brief description of the ALEPH detector is given 
is Section |21 The signal and background simulations are summarized in Section El The 
overall search strategy is presented in Section 0] and the specific selection algorithms in 
Section El The results are reported in Section El 

Table 1: Integrated luminosities and centre-of-mass (CM) energy ranges for the data collected 
by the ALEPH detector for the years 1999 and 2000. 



Year 


Luminosity (pb 1 ) 


CM energy ( GeV) 


2000 


7.3 ±0.05 
125.9 ±0.6 
81.4 ±0.4 


207 - 209 
206 - 207 
204 - 206 


1999 


42.6 ±0.2 
87.2 ±0.4 
79.9 ±0.4 
28.9 ±0.1 


201 - 203 
199.5 
195.5 
191.6 
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Figure 1: Branching fraction of benchmark fermiophobic Higgs boson (defined in the text) into 
boson pairs as calculated by HZHA [5] . 

2 The ALEPH detector 

A detailed description of the ALEPH detector can be found in Ref. and of its 

performance in Ref. ^2j- Here, only a brief description of the detector elements and 
the algorithms relevant to this analysis is given. 

The trajectories of charged particles are measured with a silicon vertex detector 
(VDET), a cylindrical drift chamber (ITC) and a large time projection chamber (TPC), 
all immersed in a 1.5 T axial magnetic field provided by a superconducting solenoidal 
coil. The energy of electrons, photons, and hadrons is measured with the electromagnetic 
(ECAL), the hadron (HCAL) and the luminosity (LCAL and SiCAL) calorimeters. The 
ECAL, placed between the TPC and the coil, is a highly segmented calorimeter, which 
is used to identify electrons and photons, and to measure their energy and position. The 
LCAL and SiCAL extend the calorimetric coverage down to 34 mrad from the beam axis. 
The HCAL consists of an instrumented iron return yoke. It provides the measurement of 
hadronic energy and, together with external chambers, muon identification. 

In the following, good tracks (or simply tracks) are defined as charged particle tracks 
reconstructed with at least four hits in the time projection chamber, originating from 
within a cylinder of 20 cm length and 2 cm radius, coaxial with the beam and centred on 
the nominal collision point, and with a polar angle with respect to the beam such that 
| cosfl| < 0.95. 

In this analysis, all searches make use of the same lepton identification criteria, where 
needed and applicable. Electrons are identified by comparing the momentum measured 
in the tracking detectors with the energy measured in the ECAL, by the depth and 
shape of the ECAL shower, and by the specific ionization information from the TPC, 
when available. Muons are identified by their characteristic hit pattern in the hadron 
calorimeter, and must have at least one associated hit in the muon chambers. Lepton 
identification is described in detail in Ref. [TBI . 
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Global event quantities, such as the total visible energy (-Etot), the total visible mass 
(M to t), or the missing energy (E m i SS ) are measured with an energy-flow algorithm [T2"] . 
which combines individual tracker and calorimeter measurements into energy-flow 
"particles". These objects are classified as photons, electrons, muons, neutral or charged 
hadrons. The jets used in the present analysis are obtained by clustering the energy- 
flow particles with the Durham jet-finding algorithm |14j . and allow various event 
topological variables to be determined. These include, for example, the acollinearity angle 
#aco (respectively the acoplanarity angle $ aco ) between two jets (respectively between their 
projections onto the plane transverse to the beam axis) in an event forced to form two 
jets, the transition values yij of the resolution parameter y cnt at which the number of jets 
in an event switches from j to i jets, or the good track multiplicity in jet i. Other 
specific variables, related either to the event topology, to the jet properties, or to the 
lepton characteristics, are described in turn in Section El 

3 Simulated samples 

Signal samples were generated using HZHA jUj for Higgs boson masses between 90 and 
117GeV/c 2 and at seven different CM energies: 191.6, 195.5, 199.5, 201.6, 204.9, 206.5 
and 208 GeV, including Higgsstrahlung and fusion processes. All decays of the Z and W 
bosons were considered. In the HZHA code, there is no spin correlation between the W 
bosons coming from the Higgs boson decay. The signal events were therefore re-weighted 
to take into account this spin correlation. The event weight was computed as the ratio 
between the full four-fermion matrix element and the HZHA matrix element [T3] . 

Event samples of all Standard Model (SM) background processes relevant for the 
Higgs boson search were also generated: the Bhabha process was simulated with 
BHWIDE 1.01 jinj, qq, dimuon and ditau events with KK 4.14 [17], 77 processes with 
PH0T02 18|. In the following, these processes will be grouped under the label "two- 
fermions" . WW production was simulated with K0RALW 1.51 ^J] and the remaining four- 
fermion processes with PYTHIA 6.1 [201 • The background event samples were generated 
at the same CM energy values as the signal. The simulated sample sizes are at least 
a factor 20 greater than the data. A detailed simulation of the detector response was 
applied to both background and signal events. 

4 Search strategy 

4.1 Event classes, topological searches and targeted channels 

The event selection is subdivided in a number of topological searches, each of which 
targets a specific final state (or channel) arising from the ZWW* production. The list of 
channels addressed in this paper is given in the second column of Table with the Z, 
the W and the W* decays given in this order, together with the corresponding branching 
fractions. Altogether, almost 80% of the possible final states are targeted by the selection 
algorithms developed for this study. 

Events are first separated in four exclusive classes according to the number and the 
energy of identified leptons in the final state. The four classes, further subdivided in 
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different topologies (or subclasses), are defined as follows, and are displayed in the first 
column of Table El 

• Fully Hadronic (Class 1): in this class, only events with neither energetic nor isolated 
identified leptons are selected. It addresses final states exclusively with hadronic jets, 
with (class la) or without (class lb) missing energy, depending on whether the Z 
decays in a pair of quarks or a pair of neutrinos. 

• Two Hard Leptons (Class 2): in this class, events with at least two energetic 
identified leptons are selected. It addresses final states in which the Z decays into a 
lepton pair. Five different topological searches were developed according to the W 
and W* decay modes. 

• One Hard Lepton (class 3): in this class, events with exactly one energetic identified 
lepton are selected. It addresses final states in which the W decays to £u. Four 
different topologies were defined according to the hadronic activity and the missing 
energy, to target the remaining Z and W* decays. 

• One Soft Lepton (class 4): in this class, events with exactly one isolated identified 
lepton that does not meet the momentum requirement of class 3 are selected, to 
address the W* — > tv decays. Two different topological searches were developed 
according to the hadronic activity and the missing energy, to address the hadronic 
and invisible Z decays. 

These definitions, and the corresponding topological searches, were developed to 
minimize the cross-channel contamination between the different subclasses. 

4.2 Topological search optimization 

In each of the topological searches, the selection criteria were tailored to optimize the 
combined sensitivity to a Higgs boson mass hypothesis of 110GeV/c 2 , which is near the 
expected experimental sensitivity in the fermiophobic scenario. To do so, the expected 
combined confidence level on the signal hypothesis, (CL S ), that would be obtained on 
average if no signal were present, is minimized [21] with respect to the position of the cuts 
on most of the selection variables. An estimate of the value of (CL S ) is determined with a 
toy Monte Carlo method using the approximate formula of Ref. [22 , with the algorithm 
of Ref. [23- 

For each topology, the determination of (CL S ) requires the expected number of 
background events (iVb), the number of signal events expected from the targeted final 
state N B , and the expected distribution of a variable D aimed at discriminating between 
the signal and the background. The observed values of the confidence levels on the signal 
and on the background hypotheses, CL S and CLt,, are obtained in the same way from the 
number of events observed (N^) and the value of D for each of these events. 

In the optimization process the values of and N s correspond to the CM energies and 
the integrated luminosity collected in the year 2000. Lower energy data sets would indeed 
not contribute significantly to the combined sensitivity to a 110GeV/c 2 Higgs boson 
signal, and are therefore absent from all distributions presented in Section El Since, 
however, they increase the sensitivity to smaller masses, the data taken in 1999 were 
included in the final result (Section EJ). 
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5 Event selection 



5.1 Preselection and class assignment 

Common preselection cuts are applied in all four classes in order to strongly reduce the 
77 and £ + £~ backgrounds. The energy within 12° of the beam axis, E%2, must be less 
than 40% of the CM energy. The acollinearity must be less than 170° for events with less 
than four tracks. Finally, the total invariant mass M tot and total transverse momentum 
Pt of the event must satisfy M tot + 6P t > 0.2>/s. 

The event-to- (sub) class assignment is based on the energies Eg. and the isolations 
I i. = 1 — cosO^t (i = 1, 2, 3) of the three most energetic leptons in the event. Here O^t is 
the angle between the zth lepton direction and the closest track in the event. If less than 
three leptons are found, the corresponding energies and isolation variables are set to and 
10 -20 respectively. Events in which the most energetic lepton is "hard" (Ee 1 > 25 GeV) 
are assigned to class 2 or 3 depending on the energy of the second most energetic lepton. 
To separate the remaining events between classes 1 and 4, a linear discriminant is 
built with Ei x , the total missing three-momentum (P m i SS ), and the isolation of the most 
energetic lepton (i^): 

D 1A = 2.3E tl + P miS s + 4.8 ln(/ £l ). (1) 

Details of the partition process are presented in Table El The criteria used to define all the 
subclasses presented in Section |4~T1 are also shown in the table. Subclasses la and lb are 
separated by a cut on the missing mass M miss . Subclasses 2a and 2t are separated from 

Table 2: Search topologies and targeted final states (with corresponding branching fractions) 
in the four event classes. For the targeted final states, the Z, W, and W* decays are given in 
this order. 



Class and topology 


Targeted Channel (BR) 


1: Fully-Hadronic 


No leptonic decay (0.422) 


la: 6 jets 

lb: 4 jets and E miss 


qqqqqq (0.328) 
i/z/qqqq (0.094) 


2: Two-Hard-Leptons 


Z leptonic decays (0.054) 


2a: plus jets 

2t: plus jets and -E m i SS 

2b: plus jets and 1 soft lepton 

2c: plus jets and 1 hard lepton 

2d: plus 1 hard lepton and 1 track 


£+£"qqqq (0.032) 
£ + £~ruqq (0.003) 
£+£~qq£u (0.010) 
£+£~£uqq (0.007) 
£+£~£u£u (0.003) 


3: One-Hard-Lepton (and E m [ SS ) 


W leptonic decays (0.171) 


3a: plus jets 

3b: plus jets and 1 soft lepton 
3c: plus 1 track and M miss 
3d: plus jets and M miss 


qq£z/qq (0.101) 
qq£u£u (0.031) 
vv£v£v (0.029) 
vv£vqq (0.008) 


4: One-Soft-Lepton 


W* leptonic decays (0.130) 


4a: plus jets 

4b: plus jets and M miss 


qqqq£z/ (0.101) 
vvqq£v (0.029) 



5 



subclasses 2b, 2c and 2d with a cut on Eg a . Subclass 2b is finally identified with a cut on 
the total hadronic energy, Eh a d- Subclasses 2c and 2d are separated from each other by 
a cut on the number of tracks (iV c h) while subclasses 2a and 2t are distinguished using 
the hadronic activity. The separation between subclasses 3a, 3b and 3c, 3d is achieved by 
cutting on the total and missing invariant mass of the event. A cut on the energy of the 
second most energetic lepton is used to distinguish subclasses 3a, 3d from 3b, 3c. Finally, 
subclasses 4a and 4b are separated from each other by cutting on the total invariant mass 
of the event. 

5.2 Class 1: Fully hadronic final state 

5.2.1 Class la: Six-jets final state 

targeted channel: ZH — > qqqqqq 

This subclass is characterized by a final state with a large number of tracks, a value of 
the total mass divided by the CM energy (M tot /y/s) close to 1, no missing longitudinal 
momentum (Pi/^/s), an intermediate sphericity (S) [21], and a high value of ln(?/34). 
Corresponding preselection cuts are applied, in addition to those applied in the four 
classes, and are presented in Table El together with the numbers of selected data, 
background and signal events. After the preselection, the background is dominated by 
WW events. 

The final selection cuts are given in the lower section of Table HJ The variables iV^ 1 ]^ 
6fi m and 6™ x are computed after having forced the event to be clustered into six jets. 
The smallest number of tracks in any jet is denoted iV™^, and the variables 8^ m and ^™ ax 
represent the smallest and largest angle between any pair of jets. The jet pair with the 
invariant mass closest to 91.2GeV/c 2 is assigned to the Z. The two least energetic jets 



Table 3: Details of the partition of events into each of the thirteen subclasses. Energies are 
expressed in GeV and masses in GeV/c 2 . 



Subclass 


Selection Criteria 






Ee 2 






-^miss 


-Shad 


iV ch 


la 


< 25 






< 13 


< 60 






lb 


< 25 






< 13 


> 60 






2a 


> 25 


> 20 


< £ 


I 




> 60 




2t 


> 25 


> 20 


< £ 






< 60 




2b 


> 25 


> 20 


> £ 


J 




> 60 




2c 


> 25 


> 20 


> £ 


J 




< 60 


> 4 


2d 


> 25 


> 20 




i 




< 60 


= 4 


3a 


> 25 


< 10 




> 0.4 


< 95 






3b 


> 25 


[10,20] 




> 0.4 


< 95 






3c 


> 25 


[10,20] 




< 0.4 


> 95 






3d 


> 25 


< 10 




< 0.4 


> 95 






4a 


< 25 






> 13 > 0.6 








4b 


< 25 






> 13 < 0.6 
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Table 4: Selection criteria for each subclass in class 1. The numbers of signal (JV B ), background 
(iVb) and data (iVd) events are given in the table for the year 2000. Energies, momenta and 
masses are expressed inGeV, GeV/c andGeV/c 2 , respectively. 



Subclass la 


Cuts 




N b 


A r d 


Preselection 


N d , > 25 










M^Vi > 0.6 










< 0.15 










0.1 < S < 0.85 










ln(j/34) > -5.3 


10.9 


1452.0 


1401 


Topological 


S > 0.13 


10.7 


1374.1 


1325 




N ch > 32 


9.17 


866.9 


845 




ln(t/ 56 ) > -7.2 


8.18 


494.5 


481 




Hvi2 + J/34 + Vm) > -0.83 


8.07 


476.4 


456 




^ > 


6.28 


231.5 


247 




(0™ in + 0!» ax ) > 190.° 


6.21 


212.3 


229 


Masses 


Mr < 117. 


6.10 


186.7 


200 




M 3 > 13. 


5.74 


163.7 


168 



Subclass lb 


Cuts 


A s 




N d 


Preselection 


cos(0 mlss ) < 0.9 
E 12 < 0.05 
N dt > 12 
AWv/i > 0.35 
P,/^<0.2 
0.03 < S < 0.8 










Hvu) > -7.6 


3.60 


190.4 


211 


Topological 


Wei, > 24 


1.77 


31.2 


46 




Mi/34) > -6.6 


1.73 


28.2 


43 


Anti-qq 


*aco < 179. 


1.69 


26.1 


38 




Pt/y/s > 0.035 


1.63 


22.0 


28 


Anti-WW 


T > 0.74 


1.33 


16.3 


23 




PiA/5<0.12 


1.25 


12.1 


19 






1.23 


11.3 


18 




ejj in > 23.° 


1.16 


8.61 


15 




(^h" + 8?r) > 170.° 


1.11 


8.18 


13 



are assigned to the W*, and the remaining two jets to the W. The invariant mass of the 
Z, W and W* are denoted Mi, M 2 and M 3 , respectively. 

The discriminant variable M2 enters the CL S computation and is shown in Figs. 
and |2b ; after the preselection and the final selection, respectively. In this subclass, the 
targeted-signal efficiency is 61% and the number of expected background events is 163.7. 

5.2.2 Class lb: Four-jets and missing transverse momentum final state 

targeted channel: ZH — > vv qq qq 

This subclass is characterized by a final state with a large missing mass. The preselection 
cuts used for subclass la are also relevant with, in general, weaker cut values, as shown 
in Table |3J The class lb preselection has two additional cuts on the value of the cosine of 
the polar angle of the missing momentum cos(6 l m i SS ) and on E\i- After the preselection, 
the proportion of main background events is 43% WW, 30% qq, 12% Wez/ and 12% ZZ. 

The final selection cuts are detailed in the lower section of Table EJ The acoplanarity 
and P t j \fs are used as well as the thrust T j21| of the event. The event is forced to form 
four jets. The smallest number of tracks in any jet is denoted A™^, and the variables 
#j™ in and #j™ ax represent the smallest and largest angle between any two jets. 

The discriminant variable M tot / \fs is shown in Figs. Eh and Eli after the preselection 
and final selection respectively. In this subclass, the targeted-signal efficiency is 32% and 
the number of background events is 8.2. 

5.3 Class 2: Final state with more than one hard lepton 

5.3.1 Class 2a: Two leptons and four jets 

targeted channel: ZH — f £ + £~ qq qq 

The dominant backgrounds after preselection are qq, semi-leptonic WW as well as ZZ 
events where one of the Z bosons decays into hadrons and the other into leptons. The 
rejection of qq and WW events is achieved by applying a cut on the variable 2/45 and 
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Figure 2: Discriminant variable M2 for class la events after the preselection (a) and after the 
final selection cuts (b). Discriminant variable M tot /^ for class lb events after the preselection 
(c) and after the final selection cuts (d). All distributions are obtained from year 2000 data. 



requiring the invariant mass of the two leptons (m^ c ) to be in a window around the 
nominal Z mass (mz). The remaining ZZ background is reduced by cutting on the angle 
between the two leptons and on their total energy. The details of the preselection and 
selection are shown in Table El 

The discriminant variable, inspired from the Higgs boson mass is then computed as: 

D x = \J (E tot - E h - E h y - (P tot - 4 - P £2 )2 , (2) 

where P4 and Pi 2 are the momenta of the leptons associated with the selected pair and P tot 
is the total measured momentum. The discriminant variable is shown in Figs. EH and Eb 
after preselection and final selection, respectively. In this subclass, the targeted-signal 
efficiency is 74%. The expected background is 0.67 events. 
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5.3.2 Class 2t: Two leptons and missing energy 

targeted channel: ZH — > rz/ qq 

This is the only case where the one- or three-prong hadronic tau decays can be 
distinguished efficiently from other hadronic decays. The rejection of qq and WW events 
is achieved by requiring the reconstructed mass of the two leptons to be in a window 
around the Z mass. A cut on the transverse momentum of the lepton pair (P t (£ + £~)) 
reduces the ZZ background. The selection cuts are shown in Table El In this subclass, 
the targeted-signal efficiency is 60%. The expected background is 0.46 events. 

5.3.3 Class 2b: Two leptons, two jets and one soft lepton 

targeted channel: ZH — > £ + £~ qq£u 

This subclass is characterized by a third lepton and significant hadronic activity. The 
selection proceeds in a similar way as in the subclass 2a by concentrating on the two 
leptons associated with the Z boson decay. The isolation for the lepton that is the most 
anti-parallel to the missing momentum (Ie A ) is used in the selection. The lepton tends 
to be more isolated in the signal than in the background. The selection cuts are shown 
in Table El After the final selection the dominant background is ZZ events decaying into 
£+£~hh. 

The pair of leptons which have the same flavour and opposite charge and give the best 
estimate of the Z mass are associated with the Z boson decay. The following discriminant 
is used: 

D 2 = y/ (E tot - E z y - (P tot - P z )2 + 2(E £a (E tot - E z ) + P £:i (P tot - P z )) , (3) 

where Pz and E% are the momenta and energy of the Z boson, determined from the two 
assigned leptons. The additional term with respect to D\ introduces the correction needed 
to take into account the undetected neutrino, assumed to be produced back-to-back to 
the third lepton. The discriminant variable is shown in Figs. Eb and Eli before and after 
the full selection, respectively. In this subclass, the targeted-signal efficiency is 71%. The 
expected background is 0.16 events. 

5.3.4 Class 2c: Three leptons and two jets 

targeted channel: ZH — > £ + £~ £v qq 

This subclass is characterized by a third lepton and low hadronic activity. Compared to 
subclass 2b, the missing transverse momentum is higher, which makes it more difficult 
to distinguish from semi-leptonic WW events. The mass window for the reconstructed Z 
mass is also broader due to the larger combinatorial background. A cut on the transverse 
momentum of the Z boson is applied to reduce the ZZ background. 

The signal and background distributions of the discriminant variable, D2, are shown 
in Figs. Eh andEF- The dominant background is ZZ events. In this subclass, the targeted- 
signal efficiency is 91%. The expected background is 0.17 events. 
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Table 5: Selection criteria for each subclass in class 2. The numbers of signal (N s ), background 
(iVb) and data (iVd) events are given in the table for the year 2000. Energies, momenta and 
masses are expressed in GeV, GeV/c and GeV/c 2 , respectively. 



Subclass 2a 


Cuts 




iVb 


N d 


Preselection 


N ch >8 


1.03 


48.8 


39 


Auti-qq, WW 


Mv*) > -7- 
m r z K — niz\ < 14. 


0.62 


6.29 


2 


Anti-ZZ 


e ( , t2 > i35.° 

E t , + E e2 < 95. 


0.54 


0.67 


2 



Subclass 2b 


Cuts 


N s 


N b 


N d 


Preselection 


N ch > 7 


0.29 


10.7 


12 


Anti-qq, WW 


ln(y 4 5) > -8. 
I m z c — m zj < 20. 
ln(/ f J > -7. 


0.19 


1.08 





Anti-ZZ 


6> Wa > 142.° 
Ei, + < 98. 


0.18 


0.16 






Subclass 2d 


Cuts 


N s 


Ab 




Preselection 


p,A/s>o.n 


0.10 


3.13 


3 


Anti-ZZ 


T < 0.98 
$ aco < 176.° 


0.09 


0.93 


1 


Auti-WW 


ln(/ t J > -9. 


0.09 


0.58 






Subclass 2t 


Cuts 


N s 


jv |3 


Ad 


Preselection 


P,/^s > 0.002 


0.08 


2.59 


2 


Selection 


l m z c — m zl < 23. 

p^+r) < 60. 


0.06 


0.46 


1 




Subclass 2c 


Cuts 


A^ 


A b 


N d 


Preselection 


P t /y/s> 0.01 

Na, > 4 


0.23 


1.73 


2 


Auti-qq, WW 


111(2/45) > -11. 
I m z c — m z| < 23. 
ln(/,J > -11. 


0.20 


0.18 





Anti-ZZ 


P tz < 60. 


0.20 


0.17 






5.3.5 Class 2d: Three leptons plus one track 

targeted channel: ZH — > £z/£z/ 

This subclass is characterized by a small branching fraction (the targeted channel 
represents 5% of the events in class 2) but with a clear topology: four leptons in the 
final state, one of them being soft. To reduce the ZZ background, cuts on the thrust and 
the event acoplanarity are applied. The remaining WW events are rejected by requiring 
that the most anti-parallel lepton with respect to the missing momentum is well isolated. 
The selection criteria are detailed in TableEl In this subclass, the targeted-signal efficiency 
is 57%. The expected background is 0.58 events. 

5.4 Class 3: Final state with one hard lepton 

In class 3, in order to reduce the 77 background, additional preselection cuts are applied 
on the total transverse momentum and on the cosine of the polar angle of the missing 
momentum cos(# miss ). Depending on the subclass, a cut on the number of tracks and/or 
on the most energetic reconstructed photon (E yi ) is required. 

5.4.1 Class 3a: One lepton and four jets 

targeted channel: ZH — > qq£i/qq 

To eliminate the 77, £ + £~ and qq backgrounds, cuts are applied to the track multiplicity, 
the angle between the hard lepton and the total momentum 9^, and the mass of the 
hard W boson, reconstructed as the invariant mass of the hard lepton and the missing 
momentum (M£ i p miss ). To suppress the remaining backgrounds, mainly WW, E^, the 
hadronic acollinearity (6^0 ( no i ept )) the thrust computed without the hard lepton (T no i ept ), 
?/45, and the total transverse momentum are used. The selection criteria together with 
the numbers of signal, background and data events are summarized in Table 13 



10 



ALEPH 



| j 4 term ions 
| | 2 fermions 
| | signal (m=1 10 GeV/c) 



50 60 70 80 90 100 110 120 130 140 

D, (GeV/c ) 

(a) 



| | 4 fermions 
| | 2 fermions 
| | signal (rr=110 GeV/c ] 



70 80 90 100 110 120 130 140 150 160 

D 2 (GeV/c 2 ) 



| j 4 fermions 
| j 2 fermions 
j | signal (rr=110GeWc 2 ] 



%0 40 50 60 70 80 90 100 110 120 130 

D 2 (GeV/c 2 ) 

(e) 



c 

9 

?, i-a 



i.e; 

1-4 r 
1.2 ^ 
1^ 
0.8^ 
0.6 - 
0.4 
0.2 




50 60 70 80 90 100 110 120 130 140 

D, (GeV/c 2 ) 

(b) 



|aia 
|ai4 

0.12 
0.1 
0.08 
0.06 
0.04 
0.02 



io 70 80 90 100 110 120 130 140 150 160 

D 2 (GeV/c 2 ) 

(d) 



& 

§0.18 
3> 

2 0.1 6 
0.14 
0.1 2 
0.1 
0.08 
0.06 
0.04 
0.02 



gr 

30 40 50 60 70 80 90 100 110 120 130 

D 2 (GeV/c 2 ) 

(f) 



Figure 3: Discriminant variable, D±, for subclass 2a after the preselection (a) and the final 
selection (b). In plots (c) and (d), the discriminant variable, D2, is shown after the preselection 
and the final selection, respectively, for events in class 2b. The two lower plots (e) and (f) show 
the same variable, after the preselection and the final selection, respectively, for events in class 
2c. All distributions are obtained from 2000 data. 



Figures. and EJd show the discriminant variable, M tot ( hadr ) , after the preselection 
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and after all cuts, respectively. In this subclass, the targeted-signal efficiency is 39%. The 
expected background is 2.9 events. 

5.4.2 Class 3b: One lepton, two jets plus one soft lepton 

targeted channel: ZH — > qq£v £v 

The selection procedure is very similar to that in subclass 3a, and is shown together with 
the numbers of signal, background and data events after each selection step in Table 
The 77, £ + £~ and qq backgrounds are rejected by cuts on the number of tracks, on 
the total acollinearity and on the transverse momentum. To further suppress qq events, 
cuts on the angle between the hard lepton and the total momentum as well as on the 
event sphericity are applied. Finally, most of the WW events are removed by a cut on 
the hard lepton energy which is sensibly smaller for signal events than for WW events, 
on the hadronic acollinearity, on y 12 ( no iept)> computed without the leptons, and on y 34 . 
Remaining background events are mainly semi-leptonic WW decays with a soft lepton 
produced in a jet. 

A good estimate of the Higgs boson mass can be obtained from twice the sum of the 
lepton energies. The following variable is therefore used as discriminant variable: 

D 3 = 2(E ei +E i2 ). (4) 

The discriminant variable, D 3 , after the preselection and after all cuts is shown in Figs.Et 
and |Ui respectively. In this subclass, the targeted-signal efficiency is 60%. The expected 
background is 0.76 events. 

5.4.3 Class 3c: One lepton and one track 

targeted channel: ZH — > vv £v £v 

This subclass is characterized by a significant acollinearity between the two leptons, and 
by a small invariant mass. Only events with exactly two tracks are kept. The 77 events are 
rejected by a cut on the acollinearity and on the transverse momentum. The total mass 
and the acoplanarity are used to reject £ + £~ events as well as part of the WW background. 
A final cut on the total missing mass and on 9^ removes most of the remaining WW and 
ZZ background. 

The discriminant variable, D 3 , after the preselection and after all cuts is given in 
Figs. and 03 respectively. The selection criteria together with the numbers of signal, 
background and data events are summarized in Table El In this subclass, the targeted- 
signal efficiency is 58%. The expected background is 3.5 events. 

5.4.4 Class 3d: One lepton and two jets 

targeted channel: ZH — ► vv Iv qq 

Events in this subclass are characterized by a single hard lepton with some soft hadronic 
activity. The 77 and £ + £~ events are rejected by cutting on the aplanarity, A [23]. The 
dominant WW background is then reduced in two stages. First, the angle between the 
hard lepton and the total momentum and the energy of both the first and second most 
energetic leptons are used. Then, a final rejection is achieved by cuts on the sphericity, 
on M tot (hadr) and on y 12 (no i ept ). 
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Figure 4: Discriminant variable, M tot (h a dr)> f° r subclass 3a events after the preselection (a) 
and after the final selection cuts (b). Discriminant variable, D3, for subclass 3b events after the 
preselection (c) and after the final selection cuts (d). All distributions are obtained from 2000 
data. 



The hadronic activity is included in the discriminant variable, defined as 

D 4 = \/(E tot + E ei y-(P tot -P £l ) 2 . (5) 

The discriminant variable after the preselection and after all cuts is shown in Figs. Eh 
andEJi respectively. The selection criteria together with the numbers of signal, background 
and data events are summarized in Table El In this subclass, the targeted-signal efficiency 
is 58%. The expected background is 0.65 events. 

5.5 Class 4: Soft-lepton final state 

5.5.1 Class 4a: One soft lepton and four jets 

targeted channel: ZH — ► qqqq£z/ 

The event selection relies on y 45 , y u r w lept ) and the lepton isolation, computed for the 
most anti-parallel lepton with respect to the missing momentum. This last variable 
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Table 6: Selection criteria for each subclass in class 3. The numbers of signal (JV B ), background 
(iVb) and data (iVd) events are given in the table for the year 2000. Energies, momenta and 
masses are expressed in GeV, GeV/c and GeV/c 2 , respectively. 



Subclass 3a 


Cuts 


N, 


N h 




Preselection 


PtA/5 > 0.05 
|cos(6» misa )| < 0.9 
K b >3 


2.7 


793.4 


823 


Anti-qq, ££, 77 


N cb > 20 
8t,s < 41.° 
l/i,P„„„ > 55. 


1.57 


37 


43 


Anti-WW 


E h < 55. 

Paco (no lept) > 137.° 
r no lept < 0.93 


1.42 


5.7 


12 


Anti-WW 


HVis) > -7-2 
Pt/s/s < 0.25 


1.24 


2.9 


5 



Subclass 3c 


Cuts 


A r s 






Preselections 


| cos(9 misa )| < 0.9 
< 15. 


0.17 


19.9 


20 


Anti-77 


A ch = 2 
9 m < 138.° 
Ptl\ft > 0.09 


0.15 


11.16 


14 


Anti-ff, WW 


M to t/Vs < 0.26 
<£„ C o < 160° 


0.14 


7.9 


9 


Anti-WW, ZZ 


M»fa > 127. 

e ts < 32.° 


0.12 


3.5 


5 



Subclass 3b 


Cuts 








Preselections 


P t /-JS > 0.05 
|cos(9 mias )| < 0.9 
N cb > 3 
E v < 40. 


0.85 


86.96 


86 


Anti-qq, M 1 77 


N cb > 5 
6U < 178.° 
Pt/^s > 0.08 


0.79 


55.40 


52 


Anti-qq 


< 41.° 
ln(5) > -2.2 


0.59 


3.63 


5 


Anti-WW 


E h < 55. 

"aco (no lept) > 135.° 
Vl% (no lept) > 0.23 

ln(i/ 3 4) > -6. 


0.49 


0.76 


1 



Subclass 3d 


Cuts 


A r a 




A r d 


Preselections 


iVV* > 0.05 

|cos(9 miss )| < 0.9 
A ch > 3 
E 71 < 15. 


1.22 


829.8 


836 


Anti-77, U 


ln(.4) > -8 


0.94 


4.92 


3 


Anti-WW 


0i,s < 57.° 
E h < 60. 
£f 2 < 6. 


0.90 


3.39 


1 


Anti-WW 


-5. < ln(S) < -2. 

Mot (hadi)/Vs < 0.17 
Mfc (no lept)) > -4. 


0.79 


0.65 






reduces the background by a factor 6 (for a 42% efficiency in signal events). A \ 2 is 
then built which takes into account the W and Z boson masses, reconstructed from the 
four jets in the event. For this, jets are paired and the mass of each pair is compared 
either to the nominal W boson mass or to the nominal Z boson mass. The jet pairing 
that minimizes the x 2 is retained: 

2 _ (M 12 - m z ) 2 (M 34 - m w ) 2 
X a 2 + ct2 W 

were M 12 and M34 are the masses of each jet pair and a is the estimated mean resolution 
on the reconstructed mass. The background is finally reduced by constraining the value of 
the x 2 i th e total hadronic mass and the lepton energy. The discriminating variable used 
for the estimation of the quoted confidence level is the reconstructed off-shell W mass, 
estimated from the missing momentum and the soft lepton momentum as follows: 

D 5 = (£w + E h ) 2 - (P miss + 4)2. (7) 

Details of the event selection are given in Table The discriminant variable, after the 
preselection and the final selection is presented in Figs. and ^p, respectively. In this 
subclass, the targeted-signal efficiency is 35%. The expected background is 4.0 events. 

5.5.2 Class 4b: One soft lepton and two jets 

targeted channel: ZH — > vv qq tv 

Events are selected by cutting on y\2 and 1/34, on the number of tracks and on the lepton 
energy. The remaining background, still 75 times larger than the signal, is reduced by 
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Figure 5: Discriminant variable, D%, for subclass 3c events after the preselection (a) and the 
final selection (b). Discriminant variable, D4, for subclass 3d events after the preselection (c) 
and the final selection (d). All distributions are obtained from 2000 data. 



additional cuts on the total hadronic mass, the lepton isolation and on the hadronic 
acollinearity. Since at least three neutrinos are expected in this channel, a cut on the 
missing mass is used in the final selection. The dominant background after the selection 
is WW. 

The discriminant variable, -D5, after the preselection and the final selection is presented 
in Figs. Efc and Eli, respectively. The expected number of events after the signal selection 
is shown in Table In this subclass, the targeted-signal efficiency is 54%. The expected 
background is 8.0 events. 



6 Results 

The numbers of signal (mu = 110 GeV/c 2 ), background and data events for the years 1999 
and 2000, for each class, are given in Table |H1 The numbers of observed events agree well 
with the expectations for all subclasses but 3a and 3c, where there is a small excess (~ 2a) 
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Table 7: Selection criteria for each subclass in class 4. The numbers of signal (JV B ), background 
(iVb) and data (iVd) events are given in the table for the year 2000. Energies, momenta and 
masses are expressed in GeV, GeV/c and GeV/c 2 , respectively. 



Subclass 4a 


Cuts 




N b 


N d 


Preselections 


A r cll > 10 

|cos(6» miss )| < 0.95 


4.6 


818 


794 


Group 1 


M»4s) > ~ 6 ' 

ln(j/34 ( no Icpt)) > -4.6 

ln(I<J > -3.5 


1.4 


7.4 


8 


Group 2 


X 2 < 0.02 

Mot (hadr)/Vs < 0.95 

£ fl > 7. 


1.17 


4.0 


2 



Subclass 4b 


Cuts 






A r d 


Preselections 


N ch > 3 
ff/Vs > 0.05 
|cos(0 mis s)l < 0.95 


2.4 


478 


485 


Group 1 


M?7l2 (no lept)) > -7- 
Mfei (no lopt)) > -2.6 

N ci > 9 
> 7. 


1.35 


99.7 


82 


Group 2 


Mt„, (hadr) 6 [50., 103.] 
ln(iij > -4.5 

^aco (no lept) > 140.° 
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Figure 6: Discriminant variable for subclass 4a events after the preselection (a) and after the 
final selection cuts (b). Discriminant variable for subclass 4b events after the preselection (c) 
and after the final selection cuts (d). All distributions are obtained from 2000 data. 
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Table 8: Numbers of signal events for a 110GeV/c 2 Higgs boson (N B ), background (iVb) and 
data (iVd) events, as well as the value of the expected and observed confidence levels for each 
subclass, for the years 1999 and 2000 together. 



Class 



la : ZH 
lb : ZH 



qq qq qq 
vv qq qq 



Class 1 combined 



2a 


ZH - 


-> £ + £ qqqq 


2t 


ZH - 


-> £ + £~~ rvqq 


2b 


ZH - 


-> £ + £~qq£v 


2c 


ZH - 


-> £ + £~ £vqq 


2d 


ZH - 


-> £ + £~ £v£v 



Class 2 combined 



3a : 


ZH 


3b : 


ZH 


3c : 


ZH 


3d : 


ZH 



qqtz/qq 
qq£v£v 
vv £v £v 
vv £v qq 



Class 3 combined 



4a : ZH 
4b : ZH 



qq qq£v 
vv qq£v 



Class 4 combined 



All combined 



6.80 ±0.06 
1.35 ±0.03 



8.15 ±0.07 



0.64 ±0.02 
0.070 ±0.006 
0.21 ±0.01 
0.24 ±0.01 
0.113 ±0.008 



1.27 ±0.03 



1.27 ± 
0.58 ± 
0.150 ± 
0.89 ± 



0.03 
0.02 
0.009 
0.02 



2.89 ±0.04 



1.41 ± 
0.85 ± 



0.03 
0.02 



2.27 ±0.03 



14.6 ±0.09 



372.7 ± 1.7 
18.0 ±0.4 



390.7 ± 1.7 



2.41 ±0.07 
1.21 ±0.08 
0.57 ±0.04 
0.31 ±0.05 
1.33 ±0.12 



5.82 ±0.17 



6.6 ±0.2 
2.01 ±0.12 

6.6 ±0.2 
0.99 ±0.15 



16.2 ±0.4 



8.92 ±0.2 
20.34 ±0.4 



29.3 ±0.5 



441.9 ± 1.9 



360 
20 



380 



5 
1 
2 

1 



9 



13 

3 

14 
1 



32 



19 



27 



448 



(CL.) 



0.60 
0.58 



0.44 



0.57 
0.94 
0.81 
0.79 
0.91 



0.43 



0.47 
0.62 
0.95 
0.44 



0.23 



0.55 
0.78 



0.50 



0.08 



CL S 



0.50 
0.79 



0.53 



0.91 
0.96 
0.83 
0.79 
0.90 



0.73 



0.84 
0.66 
0.99 
0.43 



0.47 



0.40 
0.68 



0.30 



0.26 



over the expected background (CLb=0.97). For subclass 3a, 6.6 background events remain 
after the cuts, for 1.27 signal events, while thirteen candidates are observed. Similarly, 
for subclass 3c, 0.15 signal events are expected for 6.6 background events, and fourteen 
candidates are observed in the data. The best sensitivity is achieved for subclasses with 
a low expected (CL S ). These are subclasses 3d, 3a, 4a, 2a and lb, in that order. 

The four classes are combined and the compatibility between data and background 
with and without the signal is evaluated with the log-likelihood ratio estimator In Q [23 . 
The combined expected values of the signal and background confidence levels are 0.08 and 
0.50 respectively. The observed values for the signal and background confidence levels are 
0.26 and 0.87, respectively. 

6.1 Systematic uncertainties 

The main contribution to the systematic error comes from the simulated statistics. The 
uncertainties from the production cross sections, from the simulation of the calorimeter 
energy scale and from hadronization processes are also taken into account. The beam 
background is conservatively taken into account in the analysis (via event selection 
variable ~E\%). Uncertainties from the simulated statistics are included when combining 
the searches in the several subclasses by varying the expected signal and background levels 
bin per bin in a number of toy Monte Carlo experiments, while other uncertainties are 
taken into account in a correlated way by varying the expected signal and background 
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levels coherently in the same direction for all subclasses and all bins. Each source of 
uncertainty is nevertheless handled independently, which corresponds to adding up their 
effect in quadrature. 

For the four different classes, after applying the selection cuts, the main remaining 
background comes from WW pair production. The related uncertainty, associated to 
W decays into three jets, is taken to be 2%. This is a conservative estimate from 
the uncertainty on the a s measurement The systematic errors related to the 

simulation of the absolute energy scale of the calorimeters are determined using hadronic 
Z events and are found to be ±0.9% and ±2% for the electromagnetic and hadronic 
calorimeters, respectively. The effect of a possible miscalibration of the calorimeters 
is evaluated on simulated samples by scaling the electromagnetic and hadronic part 
of the measured energy independently by these amounts. The largest of the observed 
shifts for each calorimeter is combined in quadrature. This leads to a mean uncertainty 
of the order of 3% on the remaining background level. The uncertainties originating 
from the hadronization model are evaluated by comparing WW events hadronized using 
the string model (JETSET 7.4 j2H] Monte Carlo) and with the colour dipole model 
(ARIADNE 4. 10 27 a Monte Carlo). The associated mean systematic error is 6%. Both the 
hadronization and calorimetric uncertainties are evaluated separately for each subclass. 
The uncertainties taken into account are presented in Table OH 

6.2 Combined Likelihood Ratio 

The combined log-likelihood ratio as a function of the Higgs boson mass is presented in 
Fig. El A 1.5cx excess is observed. This excess does not depend on the hypothetical mass, 
and no supporting indication of a signal is seen. An upper limit on the signal cross section 
is set as a function of the Higgs boson mass. 

Table 9: Levels of systematic uncertainties considered for each subclass. The corresponding 
total systematic uncertainty is also given in the last column, excluding the impact of the limited 
simulated statistics. 



Subclass 


Systematic uncertainty 


sources 






Calorimetry 


Hadronization 


W 


— > 3 jets 


Total 


la 


0% 


2% 




2% 


2.8% 


lb 


6% 


8% 




2% 


10.2% 


2 


0% 


0% 




2% 


2% 


3a 


2% 


4% 




2% 


4.9% 


3b 


0% 


14% 




2% 


14.1% 


3c 


2% 


8% 




2% 


8.5% 


3d 


12% 


8% 




2% 


14.6% 


4a 


2% 


0% 




2% 


2.8% 


4b 


0% 


9% 




2% 


9.2% 
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Figure 7: Log-likelihood ratio, — 21nQ, as a function of the Higgs boson mass hypothesis, mjj, 
with all data collected between 191 GeV and 209 GeV. The solid line is the result obtained from 
the data. The expected background-only and signal-plus-background likelihoods are indicated 
by the dashed lines; the light and dark shaded bands around the background expectation contain 
68% and 95% of the simulated background-only experiments, respectively. The la bands for the 
signal-plus-background hypothesis are also shown. 



6.3 Cross section upper limits 

For each class in the analysis, an upper limit on the production cross section at a given 
Higgs boson mass is derived. Figure |H1 presents the resulting 95% CL. upper limit on 
£ 2 = BR(H -> WW)a(e + e" -> Hff)/<x SM (e+e- -> Hff), as a function of the Higgs boson 
mass. The third class, providing the best compromise between a clean signature and signal 
sensitivity has the highest reach. The other three classes are nevertheless as important in 
order to get a competitive limit. Combining all four analyses, the 95% CL. upper limit on 
£ 2 as a function of the Higgs boson mass is given in Fig. El The benchmark fermiophobic 
Higgs model is drawn as the full line in the figure. The afore-mentioned excess in class 3 
results in an observed limit which is less than the expected limit. 

6.4 Mass exclusion limits 

Lower limits on the Higgs boson mass can be extracted in the context of a given model from 
upper limits on the cross section presented in the previous section. If the branching ratio 
of the Higgs boson to W bosons is taken to be 100%, the expected limit is 107.5 GeV/c 2 . 
In the fermiophobic Higgs boson scenario, it is expected that Higgs boson masses between 
97.5GeV/c 2 and 104GeV/c 2 can be excluded. For these two scenarios, due to the excess 
observed, no limit at 95% CL. can be set. 
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Figure 8: Limit on £ 2 (denned in the text) as a function of the Higgs boson mass hypothesis, 
toh, in the four different classes. The dashed line corresponds to the expected limit while the 
95% C.L. excluded region is shown by the gray area. 



However, the present search for H — > WW can be combined with the previously 
published ALEPH search for H — > 77 [3] to significantly improve the limits on the 
fermiophobic Higgs boson scenario (Fig. Q). The combined expected limit on the 
fermiophobic Higgs boson mass is 111.4 GeV/c 2 , an improvement of ~ 7GeV/c 2 on the 
sensitivity from the H — > 77 search alone [I] 1 . The observed limit is 105.8 GeV/c 2 . 

A model-independent limit can be derived by scanning the H — > 77 and H — > WW 
branching fractions. This is conveniently parametrized as: 



BRbosons 

R, 



77 



BRh^77 + BRjj— >ww + BR H 

BRh— >77/HRbosonsj 



where R 77 represents the fraction of bosonic decays into photon pairs and ranges from 
zero to one. The best limit is obtained combining the present results with those previously 



1 The actual expected limit quoted in Reference 0], 105.4 GeV/c 2 , results from using a definition of 
CL S 22 which is different from the definition adopted in the present paper. 



20 



published by ALEPH on the search for H — > 77 [I]. The 95% C.L. limit on BRbosons is 
determined at each point of the mn versus R 77 plane, resulting in the exclusion curves of 
Fig. EH 

7 Conclusions 

A search for a Higgs boson produced in association with a Z and decaying into a WW 
pair has been performed with a dataset of 453.2 pb _1 recorded by the ALEPH detector 
at centre-of-mass energies from 191 to 209 GeV. No statistically significant evidence 
for a fermiophobic Higgs boson decaying into a WW pair has been found in the data. 
Assuming Standard Model couplings to gauge bosons, a 95% C.L. upper limit on the 
ratio BR(H -> WW)(x(e + e- -> Hff)/cr SM (e + e- -> Hff) has been obtained. Combining 
this analysis with the study of 77 decays of a Higgs boson a Higgs boson mass up to 
105.8 GeV/c 2 has been excluded in the context of the benchmark fermiophobic scenario. 
A model-independent limit has been derived by scanning the H — > 77 and H — > WW 
branching fractions. This analysis complements existing searches for new physics beyond 
the Standard Model and constrains models introducing anomalous couplings in the Higgs 
sector. 
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Figure 9: Limit on £ 2 (denned in the text) as a function of the Higgs boson mass hypothesis, 
ihh. The dashed line corresponds to the expected limit while the 95% C.L. excluded region is 
shown by the gray area. The benchmark fermiophobic Higgs model branching ratio is depicted 
by the full line. 




Figure 10: The 95% C.L. limit for BRbosons as a function of mH and R 77 . The solid lines indicate 
the upper limit of exclusion regions. The crossing point between the "BRbosons = 100%" line 
and the "Fermiophobic scenario" line provides the lower limit on the Higgs boson mass in the 
benchmark scenario: tbh > 105.8 GeV/c 2 . 
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